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We analyze the results by chiral NN models for the two-nucleon system 
and calculate the predictions for the nucleon vector analyzing power of elastic 
Q . nucleon-deuteron (Nd) scattering, A y , by these models. Our conclusion is 

that a quantitative chiral two-nucleon potential does not resolve the Nd A y 
puzzle (when only two-body forces are included). 

PACS numbers: 21.30.+y, 21.45. +v, 25.10,+s, 27.10+h 

■ The term A y puzzle refers to the inability to explain the nucleon vector analyzing power 

A y in elastic nucleon-deuteron (Nd) scattering below 30 MeV laboratory energy for the 
incident nucleon by means of three-body calculations in which only two-nucleon forces are 
applied. The problem showed up as soon as it was possible to conduct three-body continuum 
calculations with realistic NN potentials. The first such calculation was performed by Stolk 
and Tjon [|IJ in 1978 using the Reid soft-core potential [0, and the first calculations with (a 
separable representation of) the Paris potential || were conducted by the Graz-Osaka group 
y—i \ in 1987 HI; in both cases, the A y predictions showed the characteristic problem. Finally, 
the 'puzzle' became proverbial when rigorous three-nucleon continuum Faddeev calculations 
using realistic forces were started on a large scale |J. Over the years, many measurements 
and calculations of Nd A y were performed (including the pd reaction that involves the 
Coulomb force ||) which all confirmed that the problem was real (see Ref. for a review): 
For energies below 20 MeV, the A y is predicted about 30% too small in the angular region 
\ around 120 deg center-of-mass angle where the maximum occurs. 

There have been many attempts to solve the problem. Already in the very early stages of 
three-body continuum calculations, when only schematic AW potentials were applied, it was 
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noticed that the Nd A y predictions depend very sensitively on the strength of the input AW 
potential in the triplet P waves |8|||] — a sensitivity that was confirmed in later calculations 
using realistic forces ||10|| . Based upon this experience, Witala and Glockle [|TT| showed in 
1991 that small changes in those 3 P wave potentials could remove the discrepancy. This 
finding gave rise to systematic investigations of the question whether the small variations of 
the low-energy phase shifts of, particularly, those triplet P waves necessary to explain the 
Nd A y are consistent with the AW data base. While Tornow and coworkers suggest 



that the low-energy AW data may leave some lattitude in the AW P waves that could 
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improve the predictions for Nd A y , Hiiber and Friar Jl3| find that it is not possible with 
reasonable changes in the NN potential to increase the Nd A y and at the same time to keep 
the two-body observables unchanged. 

Another important observation has been that conventional three-nucleon forces (when 
added to a realistic two-nucleon potential) change the predictions for Nd A y only slightly and 
do not improve them |1^,6|. Therefore, the general perception in the community has shifted 



towards the believe that the A y puzzle is the 'smoking gun' for new types of three-nucleon 



forces |15j-[18] or new physics 

However, very recently, there has been an apparent indication that the above conclu- 
sion/believe may be premature. It was reported 
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that with a two-nucleon force of a new 
type, namely, one that is based upon chiral effective field theory, the A y puzzle is resolved. 

In recent years, effective field theory methods have become increasingly popular in nu- 
clear physics. The reason for this development is the need to link conventional nuclear 
physics methods one way or the other to the underlying theory of strong interactions, QCD. 
After quark cluster models had only a limited success, it was recognized that the symmetries 
of QCD are more important than the high-energy degrees of freedom of QCD (quarks and 
gluons). The effective field theory concept distinguishes between different energy scales and 
assignes appropriate degrees of freedom for each scale while observing the over-all symme- 
tries. For traditional nuclear physics with energies below 1 GeV, the right degrees of freedom 
are nucleons and pions interacting via a force that is controlled by (broken) chiral symmetry. 
The derivation of the nuclear force from chiral effective field theory was initiated by 

l3|.[2^|. Subsequently, many 
As a result, efficient methods for 
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Weinberg |2T| and pioneered by Ordonez p2| 
researchers became interested in the subject 
deriving the nuclear force from chiral Lagrangians emerged and the quantitative nature of 
the chiral NN potential improved. This trend shows up, in particular, in the excellent 
work by Epelbaum et al. |30| where the chiral NN force was constructed using a unitary 
transformation and applying systematic power counting in next-to-leading order (NLO) and 
NNLO; and it is this potential in NLO that was applied in Ref. p0 |, resulting — seemingly — 
in a resolution of the long-standing A y puzzle. It is the purpose of this note to critically 
investigate the predictions by the chiral NN model and the implications for the Nd A y 
puzzle. 

We start our investigation by taking a close look at important phase shifts of two-nucleon 
scattering. In Table I, we list S-w&ve phase shifts and, in Table II, we show 3 P-wave phase 
shifts for energies between 1 and 200 MeV. Since charge- dependence of the NN interaction 
is not a crucial factor in the A y puzzle |12|], and since the present chiral NN potentials are 
all charge-independent and adjusted to the neutron-proton (np) data, we consider np phase 
shifts and np data. 

It is of interest to compare the phase shifts produced by the chiral NLO (next-to-leading 
order) model by Epelbaum et al. (which is the chiral potential applied in Ref. |2D| ) 
with the empirical ones from the Nijmegen multi-energy np analysis |37| (PWA93) and the 
predictions by one representative of the family of the high-precision potentials constructed 
in the 1990's (CD-Bonn |38| , |39l ). In S waves (Table I), there is, generally, good agreement 
up to 50 MeV. Above 50 MeV, differences between NLO and PWA93 show up and increase 
with energy. However, since the S'-waves are not very important for Nd A y , this may not 
have much impact on the predictions. 
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We turn now to the triplet-P waves (Table II) which are crucial for Nd A y , and focus, 
first, on the energy range below 30 MeV. The NLO 3 Po phase shifts are about 2% larger 
than PWA93, which is not significant. For 3 Pi, the differences are more drastic: The NLO 
value at 10 MeV is about 5% smaller and the one at 25 MeV is 3% smaller than PWA93. 
Finally, the NLO prediction for 3 p2 at 25 MeV is enhanced by 13%. 

To understand how variations of the 3 P phase shifts may effect observables, we consider 
the spin-orbit phase shift combination, 

A LS = ^ (-2 S3 Po - 3 Sb Pi + 5 <5 3P2 ) , (1) 

which is a measure for the strength of the spin-orbit force. Results are shown in Table III. 
At 10 MeV, one obtains A™A93 = .203 deg and = 0.212 deg for PWA93 and NLO, 

respectively, implying that the NLO value is larger by 4.4% as compared to PWA93. At 25 
MeV, the corresponding figures are Al^ 93 = 0.93 deg and Als° = 1-09 deg implying that 
NLO is 17% larger. These numbers show that the spin-orbit force of NLO is enhanced as 
compared to the Nijmegen analysis and similar enhancements are obtained when comparing 
to any of the high-precision potentials, like CD-Bonn (cf. Table III). It is well-known that 
the two-nucleon spin-orbit force is magnified in the Nd system. Therefore, a moderate 
enhancement of the NN LS'-potential leads to a substantial enlargement of the Nd A y . 
Thus, the enhanced NLO spin-orbit force as reflected in the larger values for A^s could 
very well be the explanation of the large Nd A y predictions by NLO reported in Ref. [20]. 

Next, we consider the 3 P phase shifts above 30 MeV. Here, the differences between 
PWA93 and CD-Bonn, on the one hand, and NLO, on the other, are in general larger and 
increase with energy (cf. Table II). This trend is most dramatic in 3 P2 where the discrepancies 
quickly grow into the hundreds of standard deviations. The sign of these differences are such 
as to drastically enhance the spin-orbit force, see Table III for the value at 50 MeV. Note 
that the NN t-matrix, on- and off-shell, is input to the three-body continuum calculations. 
Thus, the description of the two-nucleon data at energies above 30 MeV has impact on low- 
energy three-body predictions. Therefore, the unrealistically strong NLO spin-orbit force 
above 30 MeV may — by means of an off-shell effect — further enhance the Nd A y predictions. 
However, one may not have much confidence in this type of off-shell effect. 

The above observations trigger the question if chiral models can also make more accurate 
predictions for NN phase shifts; and if so, what are the implications for the Nd A y problem, 
in such a case. The natural way is to include higher order terms in power counting which 
should improve not only the quality of the NN phase shift reproduction but also extend the 
energy range in which it works. For that purpose, we pick up the chiral NN potential of 
Ref. |4(J (subsequently denoted by 'Idaho') that was recently developed. In the chiral Idaho 
model PHI , contact terms up to order four are included which introduces more parameters 
allowing for a better fit of the lower partial waves in a much wider energy range. In Table 
I and II, it is clearly seen that the chiral Idaho NN potential reproduces the empirical np 
phase shifts of PWA93 up to 200 MeV, accurately. 

We consider now the observable Nd A y which is the focus of this paper. We have 
calculated the predictions for Nd A y at energies 3, 10, and 65 MeV for the incident nucleon. 
The results are shown in Fig. 1. In this figure, the shaded band represents the prediction 
by the family of high-precision potentials (using always the np version of those models), 
namely, CD-Bonn jjjj, Argonne V\s ||42|| , and the Nijmegen potentials Nijm-I, Nijm-II, and 
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Nijm93 The dashed line is the predicition by the Idaho chiral NN potential |4(| and 
it is clearly seen that this prediction follows accurately the narrow band made up from the 
variations among those high-precision potentials. In conclusion, at 3 and 10 MeV, we are 
still having an A y problem if the chiral NN potential is a quantitative one. 

The evidence presented may be perceived as a convincing proof that a quantitative chiral 
potential does not resolve the Nd A y puzzle. However, there remains one objection that 
can be raised. In the literature, notably in Ref. |T2|], one can find the suggestion that the 



3 P waves at low energy are not as well determined as claimed in PWA93 [|37j. If true, 
then moderate variations of the 3 P phase shifts at low energy could be consistent with the 
low-energy NN data. This variations could be such as to enhance the low-energy spin-orbit 
force and, thus, lead to an improved prediction for Nd A y . For the purpose of seriously 
checking out this possibility, we have constructed a variation of the Idaho chiral potential 
with modified 3 P phase shifts at low energy. The column 'Modified' of Table II shows the 3 P 
phase shifts and the corresponding column of Table III reveals that for this fit the spin-orbit 
force is enhanced, similarly to NLO. However, in contrast to NLO, the 'modified' model 
is much more realistic, since the phase shifts do not diverge to unrealistic values at higher 
energies. 

We have calculated the Nd A y as predicted by the 'modified' chiral model and find, 
indeed, a considerable improvement (see dotted curve in Fig. 1). Is this the resolution 
of the A y puzzle by a chiral NN potential? To answer this question one needs to know 
if the modified chiral model is a realistic and quantitative NN potential. A precise and 
reliable answer cannot be given by just looking at phase shifts. As stressed repeatedly by 
the Nijmegen group in the past, only a direct confrontation with the NN data can reveal if a 
NN potential is quantitative or not. For that reason we have calculated the x 2 / datum for the 
fit of the world np data as represented by the 1999 database see Table IV. For a proper 



interpretation of the results of Table IV, it is necessary to establish a standard concerning 
what x 2 /datum represents a 'quantitative' reproduction of the data. This issue was debated 
a lot in the 1990 's, and the consensus that emerged was that only values below 1.1 are 
acceptable. Deliberatly, we losen this standard and consider a fit with x 2 /datum < 1.2 as 
'quantitative', while we will perceive higher values as not acceptable. 

Applying this standard, the modified model produces unacceptable values for x 2 /datum 
for all energy intervals below 75 MeV (cf. Table IV). Thus, the modified chiral model is not 
a quantitative one and, consequently, it is not the resolution of the Nd A y puzzle. Since it 
is well known that the off-shell character of the ./ViV potential plays essentially no role in 
three-nucleon scattering, one can further draw the more general conclusion: No model that 
reproduces the AW data correctly can solve that Nd A y puzzle. 

Table IV shows also the x 2 /datum of the other models discussed in this paper. It is seen 
that PWA93, Idaho, and CD-Bonn reproduce the np data below 210 MeV with the perfect 
X 2 /datum = 0.97 and 0.98. The chiral NLO potential by Epelbaum et al. |30| produces 
X 2 /datum = 37 which is grossly unacceptable. In fact, only for the interval 0-8 MeV is NLO 
acceptable. This range of validity is so tiny that no serious implications can be drawn from 
any prediction by this potential. 

Finally, we like to take this opportunity to also present an overview of other interesting 
Nd observables, which are shown in Figs. 2-8. Concerning the deuteron vector analyzing 
power, zTii, at 3 and 10 MeV (Fig. 2), it should be noted that a seemingly drastic reduction 
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of the discrepancy between theory and pd data seen at 3 MeV has its origin in large effects of 
the long range Coulomb force acting between two protons, which is not taken into account in 
our calculations 0. Taking this Coulomb force into account, iT u is also underpredicted in 
the peak region || — an equally well-known problem, which is why it would be appropriate to 
speak more generally of the vector analyzing power puzzle in elastic Nd scattering. In almost 
all cases, the Idaho chiral NN potential follows the trend of the high-precision potentials. 
The only exceptions are T20 and T21 at 65 MeV where the chiral potential predictions describe 
slightly better the data in the minimum region as compared to conventional potentials. But, 
apart from this, the quantitative chiral NN model containing contributions of higher orders 
in power counting does not produce any new signatures. 
In summary, our main conclusions are: 

• A quantitative chiral two-nucleon potential does not resolve the Nd A y puzzle on the 
two-body force level. 

• Low-energy 3 Pj NN phase shifts that "solve" the Nd A y puzzle are inconsistent with 
the low-energy NN data. 

And, finally, as a consequence of the above two points, one may expect that no quantitative 
two-nucleon force — no matter what the basis is, pure phenomenology, meson theory, chiral 
EFT, or anything — will ever solve the Nd A y puzzle. 

An accurate NN model requires to take chiral perturbation theory (%PT) to order four. 
At that order, also many three-body force (3NF) terms occur. According to the basic 
rules of xPT, all two- and many-body terms must be included for a complete calculation. 
Conventional 3NFs were shown to be ineffective for the A y problem P^,|6|. The advantage 
of xPT is that it provides a systematic scheme to generate all terms at a given order. As it 
turns out, there are several such chiral 3NF terms that were never considered in few-nucleon 



physics before fT5[-|l8| . It is natural to expect the resolution of the Nd A y puzzle from such 
new chiral 3NFs which, therefore, should be at the focus of future work in the field. 

At next-to-leading order (NLO) in xPT, there are no 3NF contributions. So, a calculation 
with a NLO two-nucleon potential and no 3NF seems to have formal validity. However, since 
at NLO the NN data can only be reproduced for Ti ab < 8 MeV, such a calculation is doomed 
to be inconclusive, from the outset, and higher order terms must be taken into account for 
any meaningful calculation. 
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TABLE I. 1 Sq np and 3 Si phase shifts (in degrees). 



'lab 



(MeV) 



PWA93 a 



NLO 6 



Idaho c 



CD-Bonn^ 



1 


62.068(30) 


62.05 


62.03 


62.09 


5 


63 63(8) 


63.85 


63.52 


63.67 


10 


59 96(11) 


60.26 


59.80 


60.01 


25 


50 90(19) 


50.92 


50.66 


50.93 


50 


40 54(28) 


39.29 


40.27 


40.45 


100 


26 78(38) 


21.68 


26.87 


26.38 


150 


16 94(41) 


7.26 


17.56 


16.32 


200 


8 94(39) 


-5.53 

3 Sl 


9.58 


8.31 


1 


147.747(10) 


147.70 


147.76 


147.75 


5 


118.178(21) 


118.29 


118.20 


118.18 


10 


102.611(35) 


102.84 


102.64 


102.62 


25 


80.63(7) 


80.69 


80.67 


80.63 


50 


62.77(10) 


61.74 


62.80 


62.73 


100 


43.23(14) 


38.65 


43.13 


43.06 


150 


30.72(14) 


21.67 


30.41 


30.47 


200 


21.22(15) 


7.20 


20.90 


20.95 


a Nijme£ 


$en multi-energy np analysis [| 


37]. Numbers in 


parentheses give 


the uncertainties in 



the last digits. 

fc Chiral next-to-leading order (NLO) potential by Epelbaum et al. [[3(J] using a Gaussian 
cutoff with cutoff mass A = 600 MeV. 

c Chiral NN potential by Entem and Machleidt ('Idaho-B' is used). 
^Reference |39| . 
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TABLE II. Triplet-P np phase shifts (in degrees). For notation see Table I. 



T lab (MeV) PWA93 a NLO Idaho Modified 6 CD-Bonn 



3 p 



1 


0.18 


0.19 


0.18 


0.18 


0.18 


5 


1.63(1) 


1.67 


1.64 


1.61 


1.61 


10 


3.65(2) 


3.72 


3.69 


3.62 


3.62 


25 


8.13(5) 


8.22 


8.19 


7.98 


8.10 


50 


10.70(9) 


10.84 


10.63 


10.28 


10.74 


100 


8.46(11) 


8.31 


8.17 


7.91 


8.57 


150 


3.69(14) 


2.52 


3.60 


3.66 


3.72 


200 


-1.44(17) 


-4.10 


-1.21 


-0.93 


-1.55 


1 


-0.11 


-0.12 


-0.11 


-0.11 


-0.11 


5 


-0.94 


-0.99 


-0.93 


-0.94 


-0.93 


10 


-2.06 


-2.16 


-2.04 


-2.09 


-2.04 


25 


-4.88(1) 


-5.03 


-4.82 


-4.98 


-4.81 


50 


-8.25(2) 


-8.32 


-8.22 


-8.57 


-8.18 


100 


-13.24(3) 


-12.66 


-13.36 


-13.86 


-13.23 


150 


-17.46(5) 


-15.94 


-17.59 


-17.85 


-17.51 


200 


-21.30(7) 


-18.86 


-21.28 


-21.18 


-21.38 


1 


0.02 


0.02 


0.02 


0.02 


0.02 


5 


0.25 


0.24 


0.26 


0.28 


0.26 


10 


0.71 


0.70 


0.72 


0.78 


0.72 


25 


2.56(1) 


2.89 


2.58 


2.78 


2.60 


50 


5.89(2) 


8.29 


5.86 


6.29 


5.93 


100 


10.94(3) 


22.61 


10.77 


11.24 


11.01 


150 


13.84(4) 


35.98 


13.72 


13.82 


13.98 


200 


15.46(5) 


44.31 


15.58 


15.33 


15.66 



"Uncertainties smaller than 0.005 degrees are not shown. 

^Modified version of the Idaho potential with enhanced spin-orbit force at low energies. 
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TABLE III. Spin-orbit phase shift combination Als (in degrees), Eq. (1), at 10, 25, and 50 
MeV for various models explained in Table I and II. 



Ti ab (MeV) 


PWA93 


NLO 




Idaho 




Modified 


CD-Bonn 


10 


0.203 


0.212 




0.195 




0.244 


0.207 


25 


0.93 


1.09 




0.92 




1.07 


0.94 


50 


2.73 


3.73 




2.73 




3.05 


2.73 


TABLE IV. 


X 2 /datum for the reproduction of the 1999 np datab 


ase [43 


| up to Ti ab 


= 210 MeV 


by various models explained in 


Table I and II. 














Bin (MeV) 


# of data 


PWA93 


NLO 


Idaho 


Modified 


CD-Bonn 


0-8 


81 


1.05 


1.05 




1.03 




1.32 


1.06 


8-17 


192 


0.86 


1.22 




0.85 




1.87 


U.91 


17-35 


292 


0.82 


2.15 




0.81 




1.27 


0.81 


35-75 


340 


1.03 


10.8 




1.09 




1.24 


1.01 


75-125 


239 


1.00 


11.6 




0.99 




1.06 


0.97 


125-183 


414 


1.06 


95.7 




1.07 




1.15 


1.03 


183-210 


141 


0.96 


111.7 




0.97 




1.18 


1.01 


0-210 


1699 


0.97 


37.0 




0.98 




1.27 


0.97 
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FIG. 1. Nucleon analyzing power A y for elastic Nd scattering at 3, 10, and 65 MeV. The shaded 
band represents the predictions by the high-precision potentials CD-Bonn [39|, Argonne V\% [42], 
and the Nijmegen potentials Nijm-I, Nijm-II, and Nijm93 [41] (using always the np versions of 
these potentials). The dashed line is the prediction by the Idaho chiral NN potential [40|. The 
dotted curve represents the result from the 'Modified' chiral potential (see text) and the dash-dot 
line is predicted by the NLO chiral potential by Epelbaum et al. 
(nd, squares), at 10 MeV from 
(pd, crosses). 



SGf. Data at 3 MeV are from |4J 
(nd, squares), and at 65 MeV from |S| (nd, squares) and [49] 
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FIG. 2. Deuteron vector analyzing power iTn for elastic Nd scattering. Energies, bands, and 
curves as in Fig. 1. The circles are pd data: at 3 MeV from [45], at 10 MeV from [47], and at 65 
MeV from |T|. 
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FIG. 3. Tensor analyzing power T20 for elastic Nd scattering. Energies, bands, and curves as in 
Fig. 1. The circles are pd data: at 3 MeV from [^], at 10 MeV from J47|J , and at 65 MeV from 



14 



0.03 



0.02 - 



o.oi - 







-0.01 



— i — 1 — r 
3MeV 



i — 1 — i — r 

/ \ 

i •••> 
/♦♦ \ 

* v 



\ 4 i 



_L 



_L 



_L 



_J I 



30 60 90 120 150 180 



0.1 



1 — i — 1 — i — 1 — r 
0.05 1 10 MeV 

o 



-0.05 - 



-0.1 



T 



T 



v 

\ • V 

J , I , I , L 



30 60 90 120 150 180 

-i 1 1 1 1 1 1 1 r- 



o - 



3 - - 1 1- 



-0.2 



l/u 

' I I 



65 MeV 41/ 

_i i i i i i i i i i i 



30 60 90 120 150 180 

6 [deg] 



FIG. 4. Tensor analyzing power T%\ for elastic Nd scattering. Energies, bands, and curves as in 
Fig. 1. The circles are pd data: at 3 MeV from [^], at 10 MeV from J47|J , and at 65 MeV from 
0. 
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FIG. 5. Tensor analyzing power T22 for elastic Nd scattering. Energies, bands, and curves as in 
Fig. 1. The circles are pd data: at 3 MeV from [^], at 10 MeV from J47|J , and at 65 MeV from 
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FIG. 6. Spin transfer coefficients K® , , and for elastic Nd scattering at 10 MeV. Bands 



and curves as in Fig. 1. The circles are pd data from [47 



17 




--L. I 



on 



-2 100 
6 



a 

-a 



XXX.S 



» lOMeV 
\ 

x 

V 

\ 

X 



\ 



It 

/ 



\ * 

\ 1 
v f 



30 60 90 120 150 180 




30 60 90 120 150 180 

6 [deg] 



FIG. 7. Differential cross section da/dQ of elastic Nd scattering. Energies, bands, and curves 
as in Fig. 1. The crosses are pd data: at 3 MeV from [50], at 10 MeV from [47], and at 65 MeV 
from 1491 1 . The circles at 65 MeV are nd data from Eql. 
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FIG. 8. Neutron-deuteron breakup cross sections for the space-star and quasi-free scattering 
(QFS) configurations at 13 MeV along the kinematical locus S. Bands and curves as in Fig. 1. 



The circles are pd data from |52|. The crosses and squares are nd data from j53| and [54], 
respectively. 
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